The role of pinning and instability in a class of non-equilibrium growth models 
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We study the dynamics of a growing crystalline facet where the growth mechanism is controlled 
by the geometry of the local curvature. A continuum model, in (2+1) dimensions, is developed 
in analogy with the Kardar-Parisi-Zhang (KPZ) model is considered for the purpose. Following 
standard coarse graining procedures, it is shown that in the large time, long distance limit, the 
continuum model predicts a curvature independent KPZ phase, thereby suppressing all explicit 
effects of curvature and local pinning in the system, in the "perturbative" limit. A direct numerical 
integration of this growth equation, in 1+1 dimensions, supports this observation below a critical 
parametric range, above which generic instabilities, in the form of isolated pillared structures lead to 
deviations from standard scaling behavior. Possibilities of controlling this instability by introducing 
statistically "irrelevant" (in the sense of renormalization groups) higher ordered nonlinearities have 
also been discussed. 
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I. INTRODUCTION 

Growth of interfaces in a strongly temperature con- 
trolled regime has recently been a subject of considerable 
technological importance [1]. The interest in this field 
has mainly been generated by the microscopic roughness 
that originates as a result of competition among different 
effects, such as surface tension, thermal diffusion and dif- 
ferent noise factors coming into play during the growth 
process. Even in equilibrium, a crystalline facet "re- 
mains practically flat until a transition temperature is 
reached, at which the roughness of the surface increases 
very rapidly" [2]. Numerous theoretical models, starting 
from Burton, et al [2] have been proposed to account for 
a detailed analysis of the height fluctuations during the 
growth process on both sides of the roughening temper- 
ature Tr [3-5]. They found a change in mobility from 
activated growth in the nucleated phase (characterized 
by low temperatures) to non-activated growth at large 
temperatures. The pinning-depinning growth model of 
Chui and Weeks [3] has later been numerically extended 
to a polynuclear growth model by Sarloos and Gilmer [6] 
confirming a growth by island formation which is contin- 
uously destroyed by any small chemical potential favor- 
ing the growth. On a different variety of non-equilibrium 
growth, phase boundaries of growing fingers emerging out 
of a competition between morphological instabilities with 
the stabilizing effect of surface tension have also been 
studied to get an idea of how systems arrive at their 
observed patterns [7]. Theoretical forays in this front 
have continued with detailed predictions on the equilib- 
rium roughening transition [8,9] which have also found 
experimental justifications in [10]. The theoretical ef- 



forts [8,9] centered around driven sine-Gordon models, 
coupled with relevant non-linearities and then looking at 
possible roughening transitions. The experimental value 
[10] of the roughening temperature compared extremely 
favorably with the theories. 

Starting from the early theoretical efforts [2,3] to 
the present day developments [1,11], numerous discrete 
[12,13] and continuum [14-16] models have been proposed 
to study both equilibrium and non-equilibrium surface 
properties. Historically speaking, most of these nonlin- 
ear models seem to fall within cither of the KPZ or the 
Lai-Das Sarma (LD) (i. e. conserved and non-conserved) 
universality class, barring a few exceptions [12,16,17]. 
The atomistic growth process of the latter type deals 
with surfaces grown by molecular beam epitaxy (MBE) 
method, whose distinguishing feature is that there the 
growth occurs under surface diffusion conditions, with 
the deposited atoms relaxing to the nearby kinks. The es- 
sential idea employed is to modify the relaxation mecha- 
nism as a locally surface minimizing curvature. To linear 
order, this mechanism is supposed to mimic the growth 
dynamics of crystalline surfaces [13,15]. 

However, the basic atomistic feature associated with 
either of these universality classes concerns a growth 
restriction on the height difference between competing 
sites. Only relatively few digressions are known, of which 
a classic example is the equilibrium restricted-curvature 
(RC) model studied by Kim and Das Sarma [18]. The 
corresponding growth rule describes a growth restricted 
on the local curvature, | V 2 /i < N, and is obeyed at 
both the growing site and its nearest neighbors where N 
is any fixed positive integer. Particles are not allowed 
to stick to a particular site if this growth restriction is 
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not met with. We adopt the non-equilibrium version of 
the growth model proposed by Kim and Das Sarma as 
the starting point of our study of the dynamics of the 
non-equilibrium growth in an MBE process. Our basic 
theoretical approach follows the standard dynamic renor- 
malization scheme employed as in [8,9]. The technical 
details are elaborated in section 3 and a brief discussion 
on the implication of these results, from the standpoint 
of non-equilibrium critical phenomena, is provided in the 
following section. 

The remaining portion of the paper is a numerical in- 
vestigation of the same problem, starting from the same 
dynamic equation as used in the analytical approach. We 
attempt a direct numerical integration of this dynamical 
equation using Eulerian discretization schemes. However, 
here we had to restrict ourselves in 1+1 dimensions due 
to the insurgence of certain generic instabilities associ- 
ated with the growth process. Here it must be stressed 
that as observed in [19] and something which has been 
strongly supported by our numerical observations too, 
the instabilities that we see, have got nothing to do with 
the discretization scheme but is, in fact, a characteristic 
feature of the growth process itself. Identical observation, 
in context to different models have also been reported in 
[20,21]. We have checked with both the deterministic and 
stochastic versions of our model and have found more or 
less identical results as in [19], although, in our case, 
contributions from a periodic pinning potential do show 
up in certain parametric limits. However, most signif- 
icantly, a continuum analysis of this same model never 
props up with any such atomistic singularities, a result, 
which again is not unlike to those found in [19]. The rea- 
sons for this is not at all difficult to understand, since, 
in any continuum approach, we arc essentially building 
upon a perturbative analysis where the tacit assumption 
is always a vanishingly small nearest-neighbor curvature 
difference (our case) or an infinitesimally small height dif- 
ference (previous studies) between neighboring sites. But 
the catch lies in the fact that the algebra of derivatives 
cannot be analytically continued to the domain of dif- 
ference operators, if the nearest-neighbor curvature dif- 
ference (or heights, depending on the definition of the 
model) is non-negligible, as was rightfully observed by 
Dasgupta, et al. This automatically implies that any 
nonlinear discrctized process might be substantially dif- 
ferent from its continuum analogue, even though it is 
derived from the same Langevin equation. However, for 
linear growth models, there is no such disagreement be- 
tween the discrete and continuum theories. 

We study all these effects in some detail and the results 
are explicitly described in section 5. In the following sec- 
tion, we discuss the technical possibilities of controlling 
these instabilities (meaning sudden growth of a hillock 
or a hump in an arbitrary site, as compared to its neigh- 
bors) in the same line as in [19-21]. In this respect, our 
emphasis is on the method of introducing "irrelevant" 
(in the sense of renormalization) higher ordered nonlin- 
ear operators to muffle these generic instabilities while 



still retaining the structure of the basic universality class 
unperturbed. Our numerical results in 1+1 dimensions 
seem to clearly indicate that the model we are consider- 
ing belongs to the Kardar-Parisi-Zhang (KPZ) [14] uni- 
versality class in the perturbative limit of the periodic 
potential. However, in the "strong coupling" limit of this 
periodic pinning potential, higher harmonics of periodic- 
ity seem to hold sway and then our perturbative analysis 
fails, although we still seem to get numerically stable val- 
ues for exponents, which certainly are quite different from 
the KPZ ones. At the present moment, we do not have 
any analytical technique to study this "strong-coupling 
limit" (in the sense of a strong periodic potential) but 
our achievement lies in the development of a stable nu- 
merical scheme. Finally we conclude with a comparative 
discussion of numerical and analytical results in section 
7. 



II. THE MODEL AND ITS VARIATIONS 

Considering two dimensional growth pertaining to a 
lattice structure where h(f, t) is the height of the inter- 
face at time t at position r, the equation goes like 

V^ = -jV i h(r,t) + ±\V 2 h(r,t) f 

2irV . .2-Kh(f,t), „ , 
-— — sin[ ±J-L]+F + R(?,t) (1) 

where rj is the inverse mobility which fixes the time scale, 
7 is the surface tension, a is the lattice constant and V 
is the strength of the pinning potential. F is a steady 
driving force with the white noise R being defined as 

< i?(r, t)R(f, t') >= 2DS 2 (f- r*)5{t - t') (2) 

From a fluctuation-dissipation theorem, in a static ver- 
sion of the model, we relate the noise-strength D to the 
equilibrium temperature, by the relation D = r)T [9,25]. 
The biharmonic term on the right hand side of the above 
equation gives the basic relaxation mechanism in opera- 
tion. The second term is the most important term for the 
restricted-curvature (RC) dynamics and has been incor- 
porated in analogy with the KPZ equation [14]. Just as 

-> 2 

the lateral term | Vh | gives the nonlinearity in a KPZ 
growth, the curvature dependence in our model is pro- 
posed to produce a leading order nonlinearity of the form 
(V 2 /i) 2 . From phenomenological considerations, since 
the discretisation scheme prohibits a high curvature, the 
constant A is negative. Although this explains the convex 
hillocks appearing on the growing surface, the restraint 
on the unbounded growth of concave structures actually 
come from the appearance of a linear curvature term in 
cqn.(7) from a renormalization of the sine-Gordon po- 
tential. This is basically the reason for which we do not 
need to start with a linear curvature term in the start- 
ing dynamical equation, unlike as in [22,23], rather it 
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is automatically generated as a product of the dynam- 
ics of the growth process [16]. Actually we can do this 
since we are not considering any step-edge barriers. With 
identical logic, since, to start with, there is no boundary 
condition on the lateral growth property of the surface, 
a KPZ-type non-linearity is not required in the system, 
although this is indeed produced after a first renormal- 
ization of the sine term. Apparently one might think that 
this makes it an " open" system in a field-theoretic sense, 
but the point of concern for us is the phenomenology of 
the growth process, where we are only concerned with 
a curvature restriction in a periodic profile, completely 
neglecting possibilities of other initial conditions. What 
the (V 2 /i) 2 term does, therefore, is to provide the first 
order non-linear fluctuations around the steadily grow- 
ing front in the RC phase. The third term comes along 
due to the lattice structures preferring integral multiples 
of h, which means that h is measured in units of the lat- 
tice constant . As to the sign of V, we see that changing 
V — > —V, as ft — > —ft, equation (1) remains unchanged 
in the absence of the curvature term. But in the presence 
of the curvature term, a negative V, simply means that 
the periodic potential will decay to zero in the long time 
limit and the system will tend towards the Govind-Guo 
universality class [27]. So, for all practical purposes, V is 
positive. F is the steady driving force required to depin 
a pinned interface while all the microscopic fluctuations 
in the growth process are assimilated in the noise term 
R. It can be shown easily that the second and third 
terms are obtained from a variant of the stationary-state 
sine- Gordon Hamiltonian 

E[h(r,t)} =JI d 2 r[l(V 2 h) 2 -Vco S (^h)} (3) 

In the following analysis of the dynamic model, we will 
try to probe the existence of a roughening transition, 
akin to that in [8,9] and have an analytical estimate 
of this transition temperature. This obviously brings 
into question the interplay of strengths between the non- 
linear term (V 2 ft) 2 and the sine-Gordon potential. Three 
subsets of our model, defined through eqn.(l), are well 
known: 

(i) A = V = F = 0; characterizing a stationary interface 
[11,12,17]. For example, in the Das Sarma-Tamborenea 
model [12] , it was found that the growth dynamics was in- 
termediate to that between the random-deposition model 
with no relaxation (RD) and the random-deposition 
model with perfect relaxation (RDR). The growth ex- 
ponent /?(<~ 0.375) in this model is intermediate between 
that of the RD (0 ~ 0.5) and the RDR (0 ~ 0.25) models. 
This sort of growth has found experimental justification 
in the works of Yang, et al and Jeffries, et al [23] . 

(ii) A = _F = 0, V =/= 0; characterizing the growth dynam- 
ics of crystalline tensionless surfaces [25]. This model in 
equilibrium depicts a roughening transition to the high 
temperature regime of the sine-Gordon model and is ex- 
pected to model the vacuum vapor deposition dynamics 
by MBE growth. 



(iii) V = F = 0, A^0;a very special case of the 
numerical simulation predicting a "local model" for den- 
dritic growth with conventional renormalization failing 
to define a stable, non-trivial fixed point [27]. However 
numerically it is possible to determine the critical expo- 
nents as has been done bby Govind and Guo too. The 
interesting open problem here is how to arrive at these 
values analytically. 

However in both the first and second situations, de- 
tailed numerical analysis have predicted a transition from 
the conserved MBE growth to the Edwards- Wilkinson 
[28] type phase characterized by the generation of a V 2 /* 
term [25,29]. The generation of this so-called "surface 
tension" term in the dynamics although surely being a 
fall-out of the renormalization of the sine-Gordon poten- 
tial in the latter model, actually demands a greater at- 
tention. Combining this idea with the proposition put 
forward by Kim and Das Sarma [18], we therefore pose 
the most general situation concerning a surface growing 
under MBE and try to ascertain the associated rough- 
ening process throughout the whole temperature range, 
with the growth essentially occurring under a surface cur- 
vature constraint. 



III. DYNAMIC RENORMALIZATION AND THE 
RESULTS OBTAINED 

In the following analysis, wc employ standard dynamic 
renormalization techniques to probe the dynamics of our 
proposed Langevin equation (1). 

We start by first integrating over the momentum shell 
A(l - dl) <| k |< A perturbatively in A and V. There- 
after we rescale the variables back in the form k — > fc' = 
(1 + dl)k, h ft' = ft and t t' = (1 - 4dl)t. The 
various coefficients follow the rescaling 77 — > 77' = 77, 7 — > 
7' = 7, A -> A' = A, V -► V = (1 + 4dl)V, F -» 
F' = (1 + 4dl)F. We set up the perturbative scheme to 
rewrite eqn.(2) in a comoving frame of reference moving 
with velocity F/rj as 

r?— /i = - 7 V 4 ft + $(h) + i? (4) 

where $(ft) = sin[^(ft + ft)] + f (V 2 ^) 2 . Now, 

following the line of Nozieres and Gallet [8[ and Rost and 
Spohn [9], we divide both h(x,t) and R as sums of two 
parts, one an average part h, R and the other a fluctu- 
ating part 5h, SR. This latter fluctuating part actually 
defines the short-wavelength components of noise and is 
defined within the momentum shell A(l — dl) < \k\ < A. 
This gives 

- 2tt ta . .2n.- Ft., 2tt 2 rj 2 . 

$ = Vsml — (ft -I (1 ^- < Sh 2 >) 

a a 77 a z 

+ \{^hf + \<(V 2 5hf> (5) 
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and 

<5$ = -^lvcos[— (h + —)}Sh + \V 2 8h.V 2 h (6) 
a z a rj 

Now averaging out the short-wavelength components 
and including the corrections (as in [30]), we finally ar- 
rive at an expression for the renormalized mode coupling 
term, 

x Goflf-f U* — 

x e -[^A 4 (t-t') + ^0(|?-r'|,t-t')] 

- *)(»■< - rOfo - rj-)] cos[^£(t - f )] 

+ [l- 2 -^[d i h{f,t)]\r i -r' i f] sin[^^(t-f)]] (7) 

where A ~ 1/a is a suitably chosen upper cut-off with 
the Green's function G(x, t) given by 

G{x,t) = J dk e^- s -n tk4 (8) 

and 

<f>(j>,x) = [ A ^[l-J (fcp)] e -^ 4 ('-O (9) 

with p = A/?, x = 7< -*^ - and A ~ 1/L, is the inverse 
of the system size. 

Turning now to the above eqn.(7), we find that start- 
ing with a structurally non-linear term in the Langcvin 
equation, the dynamic renormalization has initiated a 
mode coupling structure with quite a different compo- 
sition compared to [8,9]. The absence of the nonlinear 
(V 2 ft.) 2 term in $ SG implies that once starting with a 
restricted curvature model of growth, the lattice struc- 
ture partakes a dynamics, where starting from an initial 
curvature dominated phase, A is later renormalized to 
zero and then at large time limit, the KPZ type non- 
linearity takes over. A relative comparison of the co- 
efficients of the RC term, as given in eqn.(l) and from 
the renormalized KPZ term, obtained from eqn.(7), eas- 
ily shows that the ratio of their respective strengths is 
equal to p-a 2 , where a is the lattice constant and A and 
A' are the strengths of the RC and KPZ nonlinearities 
respectively. To begin with, at very small times, since a 
steady growth occurs in the phase, the RC term domi- 
nates. But shortly, the pinning structure takes over the 
dynamics and from there onwards, the KPZ term comes 
into the picture. Thus, although a competition ensues be- 
tween the alternate pinning and depinning forces offered 
by the (V/i) 2 and (V 2 /i) 2 nonlinearities, at short time 



scales, for all practical purposes, in the thermodynamic 
limit, the periodic nature of the potential governs and 
the resultant scaling structure shows a KPZ type behav- 
ior. In fact, the reason for which conventional RG works 
in our model, whereas an apparently simpler version of 
it, the Govind-Guo model, comes up with an unstable 
fixed point, is due to the cross-over of our model to the 
KPZ type, in the long time, large length scale limit. In 
the following section, we tackle the problem numerically, 
which reconfirms our analytical prediction. The average 
< 5h 2 > renormalizes V and correction to the one-loop 
order looks like 

dV w = J 2 ^ < (^(o)) 2 >= -iLvdi (io) 

where is the first-ordered perturbation in potential 
V and 8h^ ' is the free height fluctuation term. From an 
analysis of the above equation, we arrive at the following 
renormalized flow equation which defines the cross-over 
temperature: 

f = (4-n )C / (11) 

where U = and n = The above equation in- 

deed is a very important one and provides a quantitative 
description of the effects of curvature restriction in the 
dynamics of our model. However, since the other flow 
equations do not really have any fundamental bearing on 
our discussion, we remove them to the appendix section 
for the more technically minded reader. Eqn.(8) clearly 
shows that in the limit n — > 4, U renormalizes 
to zero which means that in the large scale, the effects 
of the lattice potential vanish and the KPZ phase takes 
over. As consistency checks, we easily see that putting 
7 = A = 0, we reproduce the results of Rost and Spohn 
[9] and fixing 7 = A = A' = 0, the Nozieres-Gallet [8] re- 
sults are obtained. The roughening transition occurs at 

Tr = -^7-, in suitable units. This sort of a roughening 
transition was also observed in an earlier work by Tang 
and Nattermann [16]. From this point, qualitatively, the 
whole dynamics follows in the line predicted by Nozicrcs 
and Gallet and Rost and Spohn, with the biharmonic and 
structural nonlinear terms muffled by the surface tension 
term and the KPZ nonlinearity respectively. This KPZ 
regime continues until the growing interface reaches the 
next crystalline facet whereon the RC dynamics again 
comes into play and the whole mechanism goes on re- 
peating itself throughout the process of non-equilibrium 
growth. Although this apparently would imply the " irrel- 
evance" of the restricted-curvature term, but the effect of 
this curvature restriction is impregnated in the value of 
the critical temperature which is twice the value obtained 
in [8,9], an outcome of the fact that n=4 in this case, as 
compared to n=2 in [8,9]. Thus the information of the 
local geometry is conveyed through the critical tempera- 
ture, and although one finds a proliferation of interaction 
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as in the present case, the RG needs not be redone. As 
an added information, we might add that this value of 
the critical temperature is not obvious from an analysis 
of the scaling dimenion of the system and one needs RG 
to arrive at this exact value. 



IV. INSTABILITY IN THE DISCRETE GROWTH 
PROCESS 



In this section we will describe in detail results 
obtained from a direct numerical integration of the 
restricted-curvature (RC) model, in 1+1 dimensions, 
and the related instabilities in the growth process itself. 
The discretization process which we employ here is a 
standard finite-difference scheme, defined in the follow- 
ing manner: 



V 2 /(r, t) = f(r + Ar, t) + f(f- Ar, t) - 2/(r, t) (12) 



and 



V 4 /(r,t) = V 2 (V 2 /(r,t) 

= V 2 [/(r + Af,t) +f(f- Ar,t) - 2/(r,i)] 

= f(r + 2Ar, t) +f(f- 2Af, t) 

- 4 f(f + Ar, i) - 4 f(r- Ar, t) + 6 /(f, t) (13) 

where all the variables h, r, t have been suitably non- 
dimensionaliscd. In solving the above equations numer- 
ically, we start from a flat surface and invoke periodic 
boundary conditions. 

The quantities of our interest are the spatial and tem- 
poral correlation functions, defined as follows: 
spatial correlation function = (\h(r + f 1 , t) — /i(r", t)\ 2 ) P ,, 
temporal correlation function = (\h(f, t + t') — h(r, t')\ ) t ,, 
averaged over f* and t' respectively. 
We calculate the usual roughness and growth exponents 
(as defined in [1]) and compare these values with the an- 
alytical prediction of the KPZ universality class, inspitc 
of starting with a curvature dependent state. 

We studied the behavior of our ID curvature depen- 
dent growth equation (eqn.(l)) using time steps of the or- 
der of At ~ 0.001 unit and lower (upto At ~ 10~ 5 ), with 
a system size of L = 1000. Later on we checked our results 
in a larger system (L = 10,000) too but the essential con- 
clusions, to be shortly elaborated, remained unchanged. 
The results were averaged over 1024 independent runs. 
We found that for small values of A = 0.2, 0.3, 0.35, the 
roughness and dynamic exponents converged exactly on 
the ID KPZ values and the results are shown in Figs. 1 
and 2. 
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FIG. 1. R.M.S. width against spatial separation in the 
log-log scale, with A = 0.2, 0.3 and 0.35. This gives a mea- 
sure of the roughness exponent a(= 0.5), which is seen to be 
exactly the KPZ value in ID, as predicted from theory. 



For the aforesaid values of A, the roughness exponents, 
ranged from 0.49 to 0.51 with a fluctuation of ±0.02 



(3=0-33 



9 O /, 0.2 
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a a /. 0.35 




Time, In t 

FIG. 2. R.M.S. width against time in the log-log scale, with 
A = 0.2, 0.3 and 0.35. This gives a measure of the growth 
exponent /3(= 0.33), which again gives the KPZ value in ID, 
as per theoretical expectations. 

and the dynamic exponents had the values between 
0.334 to 0.337 with a fluctuation of ± 0.004 among them. 
These values were obtained using V = 10 and a = 10 
as the strength of the potential and lattice spacing re- 
spectively. However, if the lattice constant is increased, 
which physically would suggest a more coarsened pic- 
ture, keeping V unchanged, the periodic potential has 
a much reduced contribution as compared to the curva- 
ture term and in this case, there is a cross-over to the 
Govind-Guo behavior [27], with a = 1.414 ± 0.002 
and (3 = 0.392 ± 0.01. Obviously, this phase persists 
with all higher values of the lattice constant, keeping V 
unchanged. With larger values of V (~ 100) or smaller 
values of 'a'(~ 1 and smaller still), the periodic potential 
holds sway and muffles the other terms. Now the higher 
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harmonics of periodicity start popping up and this dras- 
tically modifies the overall scaling structure, such that 
the KPZ scaling is totally lost. Hence, for all of our fu- 
ture discussions, we will restrict ourselves to the values 
V = 10 and a = 10 to avoid insurgence of the higher har- 
monics, which, in the process, gives a practical measure 
of the perturbative limit, with respect to the parameters 
in hand. Although, it would have been most enterprising 
to have a theory in the " strong-coupling limit" (i.e. large 
V, small 'a' and not a field-theoretic definition) too, the 
fact that analytical techniques allow us only to trod into 
the "perturbative" domain, we refrain from encroaching 
into the "forbidden" zones, even in the numerical esti- 
mates, in the limit where the dynamics is dominated by 
strong periodic potentials. 

As stated already, a numerical analysis of the system, 
in the perturbative limit, clearly shows that for small 
enough values of A, the model is in the isotropic KPZ 
universality class. This concurs with analytical predic- 
tions of the previous section too. But in this section we 
deal with one more important property of the discretized 
system, that of evidences of "instability" in the growth 
process, even in the perturbative limit. This is inacces- 
sible through analytics and also by an atomistic simu- 
lation process, since in both these cases, the possibility 
of instability is discarded by their definition themselves 
[19]. At this point, it might seem reasonable to recapit- 
ulate the idea of " singularity" as observed in the earlier 
works [19,32,33] ." Singularity" , in our case (as in the ear- 
lier references too), is defined as a sudden and a rapid 
growth of nearest neighbor height difference at some ar- 
bitrary site, eventually leading to a computer overflow. 
Although this effect was initially observed by Tu [32] but 
both he himself and Doherty, et al [33] wrongly attributed 
the effect to "numerical artifacts". In any case, in 2D, 
they analyzed the existence of a strong- coupling (field- 
theoretically) regime where the interface develops a local 
divergence in finite time, but no such thing happened 
in the conventional weak-coupling (in a field-theoretic 
sense) regime, our point of concern. Only later incisive 
forays by Dasgupta, et al [19] and by Marsili and Bray 
[20] brought out the physics in the fore. As our observa- 
tions also indicate, the effect is a generic cause, character- 
istic to most discretized growth models, certainly in our 
model too. In connection, a very interesting result has 
only recently appeared in a work by Cuerno and Moro 
[34] . From a dynamic renormalization group study, they 
have found that in general, for all dimensions, in the over- 
damped sine-Gordon model, the linear surface diffusion 
mechanism appears only as an irrelevant perturbation 
to the basic growth dynamics, although the presence of 
the sine- Gordon damping might bring about cross-overs 
characterized by specific scaling properties. This is dou- 
bly satisfying in the context of our results, in the sit- 
uation of unconserved dynamics, although without the 
presence of any curvature term. The instability basically 
arises due to the sudden development of local hillocks or 
humps, depending on the direction of the average velocity 



of growth. In our model, the determining parameters for 
this "instability" are the curvature coefficient A and the 
strength of the pinning potential V. Our investigations 
indicate that with a positive value of A, the instability as- 
sociated is of the hillock type while for a negative A, the 
instability owes its origin to the sudden development of a 
hump. Since we are looking in the "perturbative limit", 
as has been already stated, V is kept fixed to a posi- 
tive value. An idea as to the order of magnitude of the 
heights of these hills or humps is developed from a study 
of the deterministic version of our model, which predicts 
an instability if |A/i 2 - 2iKsin(2^a) - IO7/10 > 0, 
where ho is the height of the isolated hillock (or groove). 
A non-trivial solution to the above transcendental-type 
inequality indicates an inverse dependence of A and V to 
h 0l while ho increases as the lattice constant 'a' is in- 
creased. Since we are dealing with constant values of V 
and 'a' (V = a = 10), this implies that as A increases, 
ho decreases and certainly the probability of developing 
an instability goes higher accordingly. Although, no such 
analytical prediction is possible for the actual stochastic 
equation but this nature of A dependence is still retained. 
Thus starting with an initial ho satisfying the previous 
inequality, we develop an instability shortly and depend- 
ing on the sign of A, the nature (hillock or hump) of this 
instability is then decided. It should also be mentioned 
in connection, that for a flat surface, upto A = 0.39, we 
did not find the existence of any instability but for all 
values of A, starting 0.4 onwards, the generic instability 
was seen to develop. Obviously, in making this assertion, 
we again tacitly assume fixed values of V and 'a'. Here 
one might argue from the above result that if we reduce 
A sufficiently, then even in the presence of initial hillocks 
or humps at t = 0, we can exclude the possibility of in- 
stability. But one has to keep in mind the comparative 
importance of the curvature-dependent nonlinear term 
with the periodic potential and A cannot be reduced ar- 
bitrarily. It is, however, extremely significant that even 
with largest allowable values of A, the effect of curvature 
is washed out and the system converges towards the KPZ 
fixed point in the perturbative limit. As a summary of 
this section, we can convincingly argue that the insta- 
bility observed is an inherent property of the discretized 
growth process itself and cannot be removed by reduc- 
ing the integration time step At (our check was upto the 
limit At ~ 10~ 5 ), neither is this a fallout of the dis- 
cretization scheme employed. The latter claim is amply 
borne out by the fact that even if we improve on a finite- 
difference scheme by the Runge-Kutta, the results still 
remain unchanged. Before concluding this section, we 
should probably spare a word or two about the numerical 
results to expect in 2+1 dimensions. In this context, we 
have to unfortunately declare that due to the extremely 
sensitive dependence of A and V on the instability crite- 
rion, the definition of a perturbative limit (as was V = a 
= 10 in d=l+l) appeared to be rather confused. Thus 
we have no claims to any conclusive results in d = 2+1 
and as such would avoid making any substantial demands 



G 



in this regard. Thus all the numerical results are strictly 
valid for d = 1+1, although, it might not be too wishful 
to expect that the basic KPZ universality class is still 
retained in higher dimensions also, at least, so far as the 
analytical RG calculations are to be believed. 



P=0.33 



V. CONTROLLING OF INSTABILITY 



As our previous discussions have already shown that 
we are dealing with a stochastic system which has inher- 
ent properties to develop instabilities in specific paramet- 
ric limits and thereby, destroy all possibilities of observ- 
ing the statistical properties of the system at large time 
scales, it is extremely necessary to have some method of 
controlling these instabilities. In so doing, we resort to 
the resourceful references [19,21,35] dealing with noise re- 
duction techniques to suppress instabilities in the afore- 
mentioned (KPZ, Lai-Das Sarma, etc.) nonlinear dis- 
crete models. To control instability in our model, we use 
the "irrelevant (in the RG sense) operator method" as 
used in these references. The conventional way to tackle 
this is to introduce the "irrelevant variables" (in the RG 
sense) which in turn would scale off to zero, compared to 
the leading nonlincarity in the dynamical equation. Thus 
the scaling behavior of the system remains unchanged but 
in the process, these "irrelevant" operators kill the insta- 
bility lurking within the discrete model. For our case, we 
introduce an infinite series of irrelevant terms through 
the function f(|V 2 /i| ) on the right hand side of eqn.(l), 
where 



f(x) = - (1 - exp(-cx)) 



(14) 



'c' being an adjustable parameter. For x « -, f(x) 
approaches x and in the other extreme, for x » -, f(x) 
tends to a constant value 1/c. For our RC model, a value 
of c, even as small as 0.001 is sufficient to stabilize the 
system from the instabilities. However, this is the criti- 
cal limit for c, below which the function f(|V 2 /i| 2 ) loses 
its stabilizing property Incidentally, this critical value of 
c in our model is an order of magnitude lower than was 
obtained by Dasgupta, et al with more conventional KPZ 
or LD models. For any further theoretical details on this 
topic, we refer the reader to the references [19,21] and 
avoid unnecessary repetition of well-established informa- 
tion. As a graphical back-up of our claim to stabilization, 
we demonstrate Figs. 3 and 4. 



» o /.-0.6. C=0.01 
^— □ k=0.6, C=0.001 



W 

2 
DC 



Time, In t 

FIG. 3. R.M.S. width against time in the log-log scale, with 
A = 0.6 and c = 0.01 and 0.001. This gives an idea of the 
efficacy of the instability controlled growth model. 

In Fig. 3, we plot the r.m.s. width against time, for 
A = 0.6, ho = -10.0, using c values of 0.01 and 0.001. 
The linear regions in the plots preserve the expected KPZ 
scaling with j3 — 0.337. In Fig. 4, we plot the same quan- 
tities as in Fig. 3 but for A = —0.6 and ho = —10.0 with 
identical c values. 



e X=-0.6, c=0.01 
□ X=-0.6, C=0.001 




Time, In t 

FIG. 4. R.M.S. width against time in the log-log scale, 
with A = -0.6 and c = 0.01 and 0.001. Thus with both 
A and fto negative, due to the "hump"-type instability occur- 
ing here, the linear zone, towards the origin, is rather small 
with c=0.001. 



We see that again the linear regions in the plots pre- 
serve the KPZ scaling and interestingly enough, now the 
system is stabilized even with both A and ho negative, 
something which was not possible previously. However, 
as is quite evident from Fig. 4, the process of stabilization 
is not automatic in the case where both the curvature- 
controlling parameter A and ho are of the same sign. This 
happens due to the extremely strong unstable oscillations 
arising out of this combination and with c=0.001, only 
a very small linear zone (defining [3 — 0.33) is achiev- 
able. Both in figures 3 and 4, departure from the linear 
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behavior indicate the insurgence of instabilties leading to 
non-universal non-KPZ behaviors. Our numerical results 
with a larger system size, L=10,000 shows that the effect 
is even more pronounced with larger system sizes. 

We conclude this section with comments on the possi- 
bilities of multiscaling in our RC model. In all our dis- 
cussions so far, we have most deliberately refrained from 
making comments on this issue. This is because, due to 
the presence of the periodic potential term in the dynam- 
ics, which, in effect, provides a pinning mechanism to the 
growth process, there are always possibilities of appear- 
ance of higher harmonics for all sizable (meaning larger 
than the " perturbative" value) values of the potential V. 
In fact, as we have stated even, for V larger than 10 units, 
the scaling behavior is lost altogether due to the insur- 
gence of these higher ordered periodicities. Now, naively, 
this is something to be expected too, since, our analytical 
exercise strictly borders on a perturbative estimate which 
neglects possibilities of stronger coupling terms. As soon 
as this assumption is vilified, the analytical comparisons 
no longer hold good and we land on soft grounds. In any 
case, since multiscaling exists only for a very small time, 
after which it is suppressed by instabilities, we do not 
expect to lose any important dynamic information of the 
system in neglecting this particular study, in reference to 
our model. 



VI. DISCUSSIONS AND CONCLUSIONS 

Here we discuss and analyze results obtained from 
our analytical (in 2D) and numerical (in ID) stud- 
ies in the previous sections. To begin with, an ex- 
tensive dynamic renormalization group calculation of 
the restricted-curvature (RC) model in 2+1 dimensions 
proves that the system actually belongs to the isotropic 
KPZ universality class. This result has rather far reach- 
ing consequences. The result implies that in our model of 
surface growth, which is naively expected to be explicitly 
controlled by the local curvature in the pinning-depinning 
limit, in actuality, the curvature dependence is overshad- 
owed by the KPZ nonlinearity arising out of the pinning 
potential, thereby leading to a KPZ fixed point in the 
long distance, large time limit. Thus in a normal growth 
process, particles will in general have a greater tendency 
to roll-off to its nearest neighbor sites, producing a lateral 
growth, than being controlled by the convex or concave 
structures of the locality, possibly arising out of the lo- 
cal impurities. Hence we might venture to predict that 
the effect of impurities in a dynamically growing surface 
would not show its effect by a curvature dependence of 
the locality. An identical conclusion, though in a different 
situation, was arrived at by Marsili, et al [21] through a 
reparametrization argument of the basic Langevin equa- 
tion. There they convincingly argued that a contribution 
from a local metric, arising out of a curved surface, does 
not have any effect on the overall growth process itself. 



The nonlinear exercise is largely a rethinking in the line 
of Dasgupta, et al [19], obviously in context to our curva- 
ture dependent model. The first major conclusion which 
we arrive from a numerical integration of the continuum 
equation (eqn.(l)) is that, at least, in 1+1 dimensions, 
the model belongs to the KPZ universality class. The 
numerics give clear evidence of a suppression of possible 
effects of the local curvature and driving the system to a 
KPZ fixed point, exactly as was expected from the more 
generalized 2+1 dimensional model, by analytical calcu- 
lations. This means that there is an automatic switch- 
over from the curvature dependent dynamics to the lat- 
eral growth dynamics, demanded by the system itself. 
This fact has one other indirect back-up from numerical 
observations detailed in the later part of the manuscript, 
in that, at arbitrary spatial sites, a sudden increase in 
the values of the nearest-neighbor height differences was 
seen to produce instabilities and a resulting computer 
overflow. Our model also shows generic instabilities, in- 
dependent of the size of the integration time steps or the 
definition of the discretization scheme employed. How- 
ever, our model has two independent parameters defining 
the instability, which are A and V, as opposed to only a 
single one (A) in [19]. But in effect, since we are con- 
cerned only with the perturbative limit, the second pa- 
rameter in our case is a "mild" one, compared to A. An 
interesting observation in this regards is that, for all neg- 
ative V's, the system converges towards the Govind-Guo 
fixed point [27]. The reason, again, is not difficult to 
understand, since the periodic potential, with a negative 
strength, only acts as a stabilizing higher-ordered pertur- 
bation to the Govind-Guo phase and hence is rendered 
irrelevant in the long time limit. In the perturbative do- 
main, the effect of V (> 0) shows up in reducing the 
value of c m i n . As have already been stated, our numeri- 
cal value of c m i n is an order of magnitude lower than that 
is the case with more conventional KPZ or LD type of 
dynamics. 

To conclude, in discretized versions of nonlinear con- 
tinuum models we encounter generic instabilities in the 
sense that certain sites grow much faster than their 
nearest-neighbor sites, although for linear growth mod- 
els, the discrete and continuum theories tend to agree 
with each other. However, the numerical stabilization 
scheme seems to be an extremely reliable one and one 
hopes that other works in the future will bear further 
testimony to the fact. Regarding whether discretized ver- 
sions of growth models are the more fundamental ones, 
as compared to the corresponding continuum versions, 
we believe there is an enormous scope of argument in 
this matter and our stance is that of a totally undecided 
one. This is simply because although the discretized ver- 
sions of a lot of continuum equations seem to throw up 
possibilities of generic instabilities in systems which can- 
not be accessed from the coarsened dynamics, simulta- 
neously we should not forget that atomistic versions (by 
their definition themselves) do not have any instability. 
To elucidate this further, we might consider atomistic 
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models like the restricted solid-on-solid model by Kim 
and Kosterlitz [35] and the conserved version of that [36] . 
These models, belonging to the KPZ [14] and the Lai-Das 
Sarma [15] universality classes respectively do not show 
any instability, since the basic atomistic growth rules in 
these models do not allow any surface height fluctua- 
tions to exceed a prefixed finite number, eg. unity in the 
Kim-Kosterlitz model. Thus by design, these models are 
taylor-made to suppress all such instabilities that might 
arise in a growth process. What this means is that the 
coarse-grained dynamical equations might have terms ir- 
relevant in the RG sense but these terms are certainly 
not present in their atomistic analogues. The bottomline 
story that we get from this comparison of the continuum 
models with atomistic models in the same universality 
class is that any conclusions based on a " brute- force" in- 
tegration of the continuum equation is valid also for the 
atomistic case, as long as there is no instability. But 
in the unstable zone, where normal scaling breaks down, 
such a comparison is not valid anymore. Thus the issue is 
far less than settled and we would love to have evidences, 
bearing conclusive verdict on this issue, in the future. 
However, our claim of the independence of a dynamic 
growth process to any explicit curvature effects, which 
is supported by both our analytical and numerical ob- 
servations (though in different dimensions), is something 
which can be tested from simple-to-device experiments 
and we are eagerly looking forward to such experimental 
confirmation. 



dX n 



dX' _ 8tt 5 
dl 7a 5 



nA {x '\n;K) U 2 



(19) 



(20) 



IE = 4K + » X 14 nA (K) {n . K) u* (21) 
dl 47r?77 7a' 3 



dD 



_L D ^ + ^» nA (n) {n . K)u 2 (22) 
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where U = X, K = k = ^ and n = with 

A 2 1 A 2 ' aj ja 2 ' 

7' representing the coefficient of the renormalized surface 
term [31] and A' denoting the coefficient corresponding 
to the KPZ nonlinearity generated on account of renor- 
malization. A^ % > (n; k) stand as shorthand representation 
for the integrals, where i = r],X',X,K, etc. and whose 
detailed functional forms are given below: 



A^'\n;n) 



dx 



dp p 3 Jo{p) x cos 



2ir Kxp 



72 



7 



x I dp cos[— p] exp[— — —pxp A ] 



A 7 f 



xexp[— xp —2n<fr(p,x)] 



(23) 
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VIII. APPENDIX 

The complete renormalization flow equations of the 
restricted-curvtaure model, upto the second order, are 



A (v) (n; k) = ( dx ( dp p J (p) x cos[— Kx ^ ] 
Jo Jo a 7 

dp cos[— p] exp[— — — pxp 4 )] 
X A7 

xexp[-xp 2 - 2n<j)(p, x)] (24) 



A(X')f ^ f°° dx f°° >~ ~3 t • r 27i 

A {A >(n;K)= — / dp p A J Q (p) x sm[— 
Jo x Jo a 

x J dp cos[|p] exp{-~pxp 4 )} 



2?r Kxp 2 



7 



xexp[— xp —2n(j)(p,x) 



(25) 



-dl= i4 - n)U 
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dl 



70^ 



(15) 



(16) 



a(k), x r°° dx /•°° J .. T ... . .2ttKxp 2 . 

A (K) {n-n)= — / dppJ (p)xsm[ — 

Jo x Jo a 7 



x J dp cos[^-p\ exp[— — — pxp A )] 



if' 



A 7 



xexp[— xp —2n<j}{p,x)} 



(26) 



C?7 7T 4 

~dl ~ _ 6 7 o 4 



nA^\n;n) U 2 



dl 7a 4 7 v ' ; 



(17) 



(18) 



A^(n;n)= / - / dp f J 
Jo x J 



(p) X COS 



2tt Kxp 2 



a 7 



x dp cos[yp] ex p[—-T—pxp A )} 
J X A 7 

xexp[— xp 2 — 2n(f)(p,x)] 



(27) 
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Equations (8) to (15) contain the total dynamical re- 
sponse of the system around the kinetic roughening limit. 
The behaviors of these integrals, with variation in the val- 
ues of k are more or less identical to that of Fig.l in [9], 
only in our case, the tails of each of these curves are more 
sharply defined, due to the additional multiplicative ex- 
ponential factors appearing in each of our integrals. 
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